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We consider the reconstruction technique in theories with a single or multiple (phantom and/or 
canonical) scalar fields. With the help of several examples, it is demonstrated explicitly that the 
universe expansion history, unifying early-time inflation and late-time acceleration, can be realized 
in scalar-tensor gravity. This is generalized to the theory of a scalar field coupled non-minimally to 
the curvature and to a Brans-Dicke-like theory. Different examples of unification of inflation with 
cosmic acceleration, in which de Sitter, phantom, and quintessence type fields play the fundamental 
role — in different combinations — are worked out. Specifically, the frame dependence and stability 
properties of de Sitter space scalar field theory are studied. Finally, for two-scalar theories, the late- 
time acceleration and early-time inflation epochs are successfully reconstructed, in realistic situations 
in which the more and more stringent observational bounds are satisfied, using the freedom of choice 
of the scalar field potential, and of the kinetic factor. 

PACS numbers: ll.25.-w, 95.36.-fx, 98.80.-k 

I. INTRODUCTION 

The increasing amount and precision of observational data demand that theoretical cosmological models be as 
realistic as possible in their description of the evolution of our universe. The discovery of late-time cosmic acceleration 
brought into this playground a good number of dark energy models (for a recent review, see [l|) which aim at describing 
the observed accelerated expansion, which seems to have started quite recently on the redshift scale. Keeping in mind 
the possibility, which is well compatible with the observational data, that the effective equation of state parameter w 
be less than — 1 , phantom cosmological models 0, 0, [j| share a place in the list of theories capable to explain dark 
energy. Ideal fluid and scalar field quintessence/phantom models still remain among the easiest and most popular 
constructions. Nevertheless, when working with these models, one should bear in mind that such theories are at best 
effective descriptions of the early/late universe, owing to a number of well-known problems. 

Even in such a situation, scalar field models still remain quite popular candidates for dark energy. An additional 
problem with these theories-which traditionally has not been discussed in depth-is that a good mathematical theory 
must not be limited to the description of a single side of the cosmic evolution: it should rather provide a unified 
description of the whole expansion history of the universe, from the inflationary epoch to the onset of cosmic acceler- 
ation, and beyond. Note that a similar drawback is also typical of inflationary models, most of which have problems 
with ending inflation and also fail to describe realistically the late-time universe. 

The purpose of this work is to show that, given a certain scale factor (or Hubble parameter) for the universe 
expansion history, one can in fact reconstruct it from a specific scalar field theory. Using multiple scalars, the 
reconstruction becomes easier due to the extra freedom brought by the arbitrariness in the scalar field potentials and 
kinetic factors. However, there are subtleties in these cases that can be used advantageously, and this makes the study 
of those models even more interesting. 

Specifically, in this work we overview the reconstruction technique for scalar theories with one, two, and an arbitrary 
number, n, of fields. After that, many explicit examples arc presented in which a unified, continuous description of 
the inflationary era and of the late-time cosmic acceleration epoch is obtained in a rather simple and natural way. 

The organization of this paper is as follows: in Sec. II we consider a universe filled with matter and show that it 
is possible to obtain both inflation and accelerated expansion at late times by using a single scalar field. Realistic 
examples are worked out in order to illustrate this fact. Sec. Ill is devoted to the theory of a scalar coupled non- 
minimally to gravity through the Ricci curvature. Late-time acceleration is explicitly discussed in this model, using 
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two specific examples. In Sec. IV we discuss the issue of reconstruction for a non-minimally coupled scalar field theory, 
including Brans-Dicke theory, again clarifying this case with the help of an example. Sec. V deals with de Sitter space 
in scalar-tensor theory and studies issues relevant for the conformal transformation to the Einstein frame, arriving 
to stability conditions that are important in the study of the future evolution of the late-time accelerated era. In 
Sec. V, we discuss the case of several scalar fields, beginning with the case of two scalar fields minimally coupled to 
gravity (such models have been used, e.g., in reheating scenarios after inflation). We study an explicit example and 
then the general case of n scalar fields. Sec. VII addresses our final goal, namely the reconstruction of inflation and 
cosmic acceleration from a scalar field theory, by means of a two-scalar model that reproduces cosmological constraints 
in each epoch. Further, the cosmic acceleration is reproduced with a pair of scalar fields plus an ordinary matter 
term, where the observed cosmological density parameter (Ode — 0-7) and the equation of state (EoS) parameter 
(w DE ~ — 1) are actually reproduced. Finally, Sec. VIII contains the conclusions. 



II. UNIFIED INFLATION AND LATE TIME ACCELERATION IN SCALAR THEORY 



Let us consider a universe filled with matter with equation of state p m — w m p m (here w m is a constant) and a 
scalar field which only depends on time. We will show that it is possible to obtain both inflation and accelerated 
expansion at late times by using a single scalar field <f> (see also [f| and Q). In this case, the action is 



S 



dx \/—g 



-^-R - -uifid^t - V(<P) + L, n 



(1) 



where k 2 = 8ttG, V(<fi) being the scalar potential and ui(<p) the kinetic function, respectively, while L m is the matter 
Lagrangian density. Note that for convinience the kinetic factor is introduced. At the final step of calculations, 
scalar field maybe always redefined so that kinetic factor is absorbed in its definition. As we work in a spatially flat 
Friedmann-Robertson- Walker (FRW) spacetime, the metric is given by 



ds 2 = -dt 2 + a 2 {t) ^2 dx l 



The corresponding FRW equations are written as 

^2 



{pm + P<p) > 



H 



with and p§ given by 



Pt = \u(ft>) 2 +F(0) 



{pm + Pm + P4> +P4>) 



= -wfo) 2 - V{cf>) . 



Combining Eqs. ([3]) and ([4|, one obtains 
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U>{4>) 4> 2 = -H - (p m + p m ) 



V{<f>) = 



3ir 



H - 



As the matter is not coupled to the scalar field, by using energy conservation one has 

p m + 3H{p m + p m ) = , p<j, + 3H (p^ + pj,) = . 
From the first equation, we get p m — p m oa~ 3< - 1+Wm \ We now consider the theory in which V{(f>) and w(</>) are 



(2) 

(3) 

(4) 

(5) 
(6) 



K 



l r 



3/(0) 2 + /'(</») 



: Foe -3(i+» m )F(0) j (7) 



where f{4>) = F'(cf>), F is an arbitrary (but twice differentiable) function of 0, and Fq is an integration constant. 
Then, the following solution is found (see [f| 0, H, 0] ) : 



which leads to 



<P = t, H(t) = f{t) 



a{t) = a e { >, a = — - 



3(l + ™ m ) 



(8) 



(9) 
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We can study this system by analyzing the effective EoS parameter which, using the FRW equations, is defined as 

p 2H 

- cff = - = -i-^, (10) 

where 



P = Pm + P<f> i P= Vrn + P<j> • 

Using the formulation above, one can present explicit examples of reconstruction as follows. 



(11) 



A. Example 1 

As a first example, we consider the following model: 

m = hi U 



i 



i 



(12) 



Using the solution ©, the Hubble parameter and the scale factor are given by 

t + t 



1 1 



a(t) = a Q 



tn-t 



■ arctan ~ 



(13) 



As one can see, the scale factor vanishes at t = —to, so we can fix that point as corresponding to the creation of the 
universe. On the other hand, the kinetic function and the scalar potential are given by Eqs. ([7|), hence 



u>{cj>) 



V{<f>) 



8 K{t\ + ti)U 



t 2 i+t 2 



k 2 {t\ + ^Y{tl-^Y 
hl{t\ + tl) 



(w m + l)F e- 3( ^ +1)fW 



where Fq is an integration constant and 



3^ + ^)+40U- 



j.2 _i_ j.2 



h l arctan * 
ti ti 



—F Q e-sM^W 



Then, using Eq. (JTDJ) , the effective EoS parameter is written as 

8 t{t-t + ){t + t-) 



Wcff 



-1 - 



3h 2 (t\ + tlY 



(14) 



(15) 



(16) 



where t± = ±y ° 2 1 . There are two phantom phases that occur when t_ < t < and t > i + , and another two 
non-phantom phases for —to < t < t_ and < t < i+, during which w c s > — 1 (matter/radiation-dominated epochs). 
The first phantom phase can be interpreted as an inflationary epoch, and the second one as corresponding to the 
current accelerated expansion, which will end in a Big Rip singularity when t — to. Note that superacceleration (i.e., 
H > 0) is due to the negative sign of the kinetic function u>{4>), as for "ordinary" phantom fields (to which one could 
reduce by redefining the scalar <f). 



B. Example 2 

As a second example, we consider the choice 
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We take Ho and Hi to be constants and t s as the Rip time, as specified below. Using (J7J, we find that the kinetic 
function and the scalar potential are 



w(0) = --5 



2i?i 



1 



(w m + 1)F {t s - (f>) 



3(l+w m )H 



1) 



Hi ( H 



- 2 



- 1 



^0 (t s - 4>) 



exp 

3(1 + 



3(l + w m )Hi 



■)go e 3 ' 1+ 7' H - 



respectively. Then, through the solution |J9]), we obtain the Hubble parameter and the scale factor 



H(t) 



t. 



t + ^ , a(t) = a (t, - t)- ff ° e"^ 



(18) 



(19) 



Since a(t) — > + for t — > 0, we can fix £ = as the beginning of the universe. On the other hand, at t = t s the universe 
reaches a Big Rip singularity, thus we keep t < t s . In order to study the different stages that our model will pass 
through, we calculate the acceleration parameter and the first derivative of the Hubble parameter. They are 



H 



Ho 



(t s - ty 



2Hx 



H 2 +H 



H 



(t s -ty 



;(H + 1) 



Hi ( Hi IHi 2i^o 



t 2 



t 2 



t 



ts-t 



(20) 



As we can observe, for t close to zero, a/ a > 0, so that the universe is accelerated during some time. Although this 
is not a phantom epoch, since H < 0, such stage can be interpreted as corresponding to the beginning of inflation. 
For t > 1/2 but t <C t s , the universe is in a decelerated epoch (a/a < 0). Finally, for t close to t s , it turns out that 
H > 0, and then the universe is superaccelerated, such acceleration being of phantom nature and ending in a Big Rip 
singularity at t = t s . 



C. Example 3 



Our third example also exhibits unified inflation and late time acceleration, but in this case we avoid phantom 
phases and, therefore, Big Rip singularities. We consider the following model: 



Hi 



(21) 



where Hq and H\ > are constants and n is a positive integer (also constant). The case n = 1 yields an initially 
decelerated universe and a late time acceleration phase. We concentrate on cases corresponding to n > 1 which gives, 
in general, three epochs: one of early acceleration (interpreted as inflation), a second decelerated phase and, finally, 
accelerated expansion at late times. In this model, the scalar potential and the kinetic parameter are given, upon use 
of Eqs. ([7]) and ([2T|). by 



u(<f>) 



V{4>) = 



2 nHi 



K 2 (f) n 

1 3 



|-(^ + l)F e- 3 ^ +1) (^-^-), 



(£T o n/2 +#i) nH x 



kn-l 



Then, the Hubble parameter given by the solution ([9]) can be written as 



(22) 
(23) 



H 

H(t) = H + — , a(t) = a exp 



Hot- 



Hi 



(n - l)t n 



(24) 



We can fix t = as the beginning of the universe because at this point a — > 0, so t > 0. The effective EoS 
parameter (fT0|) is 



2nH 1 t n ~ 1 
(H t" + H x ) 



(25) 
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Thus, when t — > then w e s — > — 1 and we have an acceleration epoch, while for t — > oo, iu c ff — > — 1 which can be 
interpreted as late time acceleration. To find the phases of acceleration and deceleration for t > 0, we study a/a, 
given by: 



H + H 



nHi 



H 



Hi 



(26) 



For sufficiently large values of n we can find two positive zeros of this function, which means two corresponding phase 
transitions. They happen, approximately, at 



1±*/1 



4H 



2/n 



(27) 



so that, for < t < the universe is in an accelerated phase interpreted as an inflationary epoch; for t_ < t < t + 
it is in a decelerated phase (matter/radiation dominated); and, finally, for t > t+ one obtains late time acceleration, 
which is in agreement with the current cosmic expansion. 

We now consider how the exit from inflation could be realized. First, we should note that if one refines the scalar 
field as if = dcj>y/ w(</>) in the non-phantom phase or ip = J d<j)y/—ui(<p) in the phantom phase, then the action (1) 
has the following form 



S = / d A xJ~i 



±R T ±dr<pd>><p-v(sp) 



Sn 



(28) 



Here V(<p) = V (<fi(<p)). The minus sign corresponds to the non-phantom phase and the plus sign to the phantom phase. 
For example, in the model (21) when <p is small (and therefore in the early universe), we find ip — —1\J1nH.\j (n — 
l)^^™ -1 )/ 2 and V(ip) cx ^ 4 + 4 / !) _ We should note that in the early universe (0 — » 0), <p is large. Then, the potential 
V(ip) is of the slow roll type. When the value of ip is large, the energy of the vacuum is large, which generates inflation. 
The value of ip slowly becomes small and the energy of the vacuum, and with it the curvature, decrease. When the 
curvature becomes small enough inflation could stop. For the successful exit from inflation we may need one more 
scalar field, as in the hybrid inflation scenario [irj| |. In the model (21), by construction, inflation ends as a purely 
classical theory when a = a(H 2 + H) = 0, that is, for t = t e satisfying 



: 



Hi 



nHi 



(29) 



If we include, however, quantum effects the value of the scalar field ip jumps to a larger value — as a consequence of the 
quantum fluctuations — and the vacuum acquires higher energy, what generates inflation again. In order to suppress 
the probability of those effects, in the hybrid inflation model at least one more scalar field and its coupling with ip 
had to be introduced. 



D. Example 4 

As our last example, we consider another model unifying early universe inflation and the accelerating expansion of 
the present universe. We may choose /(</>) as 



Hi + Hice 2 ^ 
1 + ce 2 



m = v::z, > ( 3 °) 



which gives the Hubble parameter 



m . BJ^l . (31) 

1 + ce ZQt 



Here Hi, Hi, c, and a are positive constants. In the early universe (t — > — oo), we find that H becomes a constant 
H — > Hi and at late times (t — > +oo), H becomes a constant again H — > Hi. Then Hi could be regarded as the 
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effective cosmological constant driving inflation, while Hi could be a small effective constant generating the late 
acceleration. Then, we should assume Hi Hi. Hence, if we consider the model with action 



S = I d xV=5 < — -r 



R uj{4>) 



>{<£) 



v{4>) 



d^d^ - v{<t>) 



f W = 2a(H i -Hi)ce 2 ^ 
k 2 k 2 (1 + ce 2a <t>f ' 

3/(<?f>) 2 + f'{4>) 3H 2 + {GH.Hi - 2a {H, - Hi)} cc 2a <>> + c 2 H 2 e ic 



™ k 2 (1 + ce 2a <t>) 

we can realize the Hubble rate given by (T3T|) with <j> = t. If we redefine the scalar field as 



J k V a 



the action S in (|32p can be rewritten in the canonical form 

S-- 

where 

H 2(H t -H l ) 



V(<p) = V(4>) = 



K 2 1 



J , n 2 



(32) 



(33) 



(34) 



(35) 



One should note that (j> — > — oo corresponds to — > and V ~ 3H 2 /k 2 , while — ► +oo corresponds to — > oo and 
V ~ 3H 2 /k 2 , as expected. At early times (<£ — > 0), V^(</?) behaves as 



3i? 2 f K 2 a(3iJ 4 + a) 



I - 



3H 2 



v 2 + o(^ 2 ) 



(36) 



At early times, <f> < and therefore, from Eq. (f3"3"| . we find K(pWa/2 (Hi — Hi) c <C 1, from which it follows that 

1 V'Up) 4a 2 /c 2 (3Hi + a) 2 , 8a (3Hi + a) 2 (H t - Hi) c 

-<p < 



3k 2 V{<pf 



3k 2 V(<p) 



27Hf 

2 (3H t + a)a 
9H 2 



27Hf 



(37) 



Then, if a <C Hi, the slow-roll conditions can be satisfied. 

We may include matter with constant EoS parameter w m . Then u>(<fi) and V(4>) are modified as 



V(ff>) 

g{4>) 



fW _ W ™ + 1 3( i +u)m)fi 
,2 o "0 



2 

K 2 



o 1 ^ 9o e 



#/(</») = + g ' - gi In (c + e" 2 ^ 



The matter energy density is then given by 

Pm = p a-^ 1+w ^ = goe-W^iM = g (c + e -2o^-3(i+- m )(^-^)/2 e _ 3 (i +U)m )^ t 



In the early universe t — ■> — oo, p m behaves as 



Pm ~ goe 



3(l+w m )(2a-Hi)t 



(38) 

(39) 
(40) 
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On the other hand, the energy density of the scalar field behaves as 

P, = ^ 2 + V^)-^ 2 + V W ^ 3 -^. (41) 

Then, if 2a < Hi (and w m > —1), the matter contribution could be neglected in comparison with the scalar field 
contribution. 

Now let the present time be t = t$. Then, we find that 

2 

Q _ K Pm 



m 2 



(42) 



and = 1 — fl m . If we assume ato 1, we find 



..2, 



ft m ^ K 90 c -3(l +Wm )(H t -H l )/2a c -3(l+w m )H,t a _ /43s 

Hence, we may choose the parameters so that J7 m ~ 0.27, which could be consistent with the observed data. This 
model provides a quite realistic picture of the unification of the inflation with the present cosmic speed-up. 



III. ACCELERATED EXPANSION IN THE NON-MINIMALLY CURVATURE-COUPLED SCALAR 

THEORY 



In the preceding section we have considered an action, JT]), in which the scalar field is minimally coupled to gravity. 
In the present section, the scalar field couples to gravity through the Ricci scalar (see 11] for a review on cosmological 
applications) . We begin from the action 



S = I d i xy / —g 



(44) 



where f(<j>) is an arbitrary function of the scalar field (p. Then, the effective gravitational coupling depends on cj>, as 
K e ff = k[1 + /(<^)] -1 ' 2 . One can work in the Einstein frame, by performing the scale transformation 

9liV = [1 + /(<£)]- % v . (45) 
The tilde over g denotes an Einstein frame quantity. Thus, the action (|4"4")l in such a frame assumes the form [l2| 



S = 



R 
2^2 



2(1 + /(</>)) K 2 (l + f(<j>)) 



v{4>) 



(46) 



The kinetic function can 



be written as W(<p) — + R s (i+f ( ) > an< ^ the ex tra term in the scalar potential 



can be absorbed by defining the new potential U(<f>) 
frame, namely 



[i+/(0)] 2 



so that we recover the action ([T]) in the Einstein 



S = I dx 4 ^ [ 4 - \ W{ct>) d^d^ - U{cf) 



(47) 



We assume that the metric is FRW and spatially flat in this frame 



dt 2 = -dP + a 2 (t)J2 dx 



i 

i ) 



(48) 
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then, the equations of motion in this frame are given by 

.2 



H 



d 2 ^ 
dt 2 



K 

-3# 



(P4>+P<t>) > 
1 



WW 



2tT'(0) 



0, 



2VF(<?!)) " vr/ V tft 

where p0 = \W{4>)<p 2 + U(<f>), p = \W{(j>)^ 2 - U(<f>), and the Hubble parameter is H = §||. Then, 

VK(0)0 2 = -45 , £/(</>) = 6H 2 + 2H. 



(49) 
(50) 

(51) 
(52) 



Note that H > is equivalent to W < 0; superacceleration is due to the "wrong" (negative) sign of the kinetic energy, 
which is the distinctive feature of a phantom field. The scalar field could be redefined to eliminate the factor W(4>), 
but this would not correct the sign of the kinetic energy. 
If we choose W((j)) and U(dj) as uj{4>) and V((j)) in ([7]), 



1 



3 ff (0) 2 + «/(</>) 



W{4>) = --g\4>) , U{4>)- , 
by using a function g(4>) instead of f(4>) in Q, we find a solution as in ([5]), 

<t> = t, H(t)=g(t). 

In (|53[) and hereafter in this section, we have dropped the matter contribution for simplicity. 
We consider the de Sitter solution in this frame, 



H = Ha = const. 



a{t) = age 



H t 



(53) 



(54) 



(55) 



We will see below that accelerated expansion can be obtained in the original frame corresponding to the Einstein 
frame (|47| with the solution (|55ll . by choosing an appropriate function /(</>). From (|55|) and the definition of W((j>) 
and U((j)), we have 



W{(j)) = -> u(<t>) 



[1 + f (</>)] K 2 



df(<f>) 



U{cj>) 



G 



2 H 



V{ct>) 



G 



(56) 



Thus, the scalar field has a non-canonical kinetic term in the original frame, while in the Einstein frame the latter 
can be positive, depending on W ((/)). The correspondence between conformal frames can be made explicit through 
the conformal transformation (|45|) . Assuming a spatially flat FRW metric in the original frame, 



ds 2 = -dt 2 + a 2 (t) ^2 dx i ) 



then, the relation between the time coordinate and the scale parameter in these frames is given by 

dt 



t = 



l (t) = [l + f(t)]- 1 / 2 a(t) 



([1 + /W] 1/2 

Now let us discuss the late-time acceleration in the model under discussion. 



(57) 



(58) 



A. Example 1 

As a first example, we consider the coupling function between the scalar field and the Ricci scalar 

1 — a<b 



a(f) 



(59) 
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where a is a constant. Then, from (|56|) . the kinetic function cu((f>) and the potential V{<j>) are 

w(</>) = 7a , = 



K 2 a 2 3 

respectively. The solution for the current example is found to be 



(60) 



, . . ~ 1 / 3a 



2/3 



_ /3a 



1/3 



exp 



a \ 2 



2/3' 



(61) 



We now calculate the acceleration parameter to study the behavior of the scalar parameter in the original frame, 



2 1 ~ / 2 
'9^ + M^ 



1/3 



2\ 1/3 1 



(62) 



We observe that for small values of t the acceleration is negative; after that we get accelerated expansion for large t; 
finally, the universe ends with zero acceleration as t — > oo. Thus, late time accelerated expansion is reproduced by 
the action (|4"4")l with the function /(</>) given by Eq. (|55|) . 



B. Example 2 

As a second example, consider the function 

/(«£) = 0-i o • 

From (|58p . the kinetic term and the scalar potential are, in this case, 

3 1 



. 2 (1 + 0- M 
The solution in the original (Jordan) frame reads 



v{4>) = ^-{i + 4>-t ) 



2a 



= T + *o - 1 , = — exp 



flo I j + *o - 1 



and the corresponding acceleration is 



a _ 1 
a ~ t 2 " 



-1 +1 



(63) 



(64) 



(65) 



(66) 



Notice that this solution describes acceleration at every time t and, for t — ► oo, the acceleration tends to a constant 
value, as in de Sitter spacetime, hence similar to what happens in the Einstein frame. Thus, we have proved here that 
it is possible to reproduce accelerated expansion in both frames, by choosing a convenient function for the coupling 

IV. RECONSTRUCTION OF NON-MINIMALLY COUPLED SCALAR FIELD THEORY 



We now consider the reconstruction problem in the non-minimally coupled scalar field theory, or the Brans-Dicke 
theory. We begin with the same scalar-tensor theory with constant parameters 4>o and Vo '■ 



S = I d A 



V{$) = Vo e" 2 ^ 



which admits the exact solution 



= 00 In 



t 

h 



H = 



2t 



2V 



(67) 



(68) 
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We choose 4>qK > 2/3 and Vq > so that t\ > 0. For this solution, the metric is given by 



ds" = —dt + On I — 



2 ,2 o 

£ \ « <Po 13 



which can be transformed into the conformal form: 



7 2 2 

as = a n 



TO 



Here 



T 
TO 



EM > 



-dr 2 + (dx i 

i=i 



TO 



( ^4>l ,\ 



and therefore, by using (|68p . one finds 



r 
to 



to 



e 



We now consider an arbitrary cosmology given by the metric 



dS"=f(r)^-dr 2 +j2(dxf^ , 



where r is the conformal time. Since 
1 



ds' 



f(r) u W 



051 — ) -V 1 - 1 ds 2 = e^ds 2 , 



/ -0 £ 



e 2 



if we begin with the action in which in (|67[) is replaced by e ip g^ Vl 

2" 



5 = / d 4 xy^e' p( 
we obtain the solution (|75|) . 



R 1 



1 - 3 



dp 



A. Example 

By using the conformal time r, the metric of de Sitter space 

ds * = _ dt 2 + e 2H t (da 

can be rewritten as 

ds 2 = 



(69) 



(70) 



(71) 



(72) 



(73) 



(74) 



(75) 



Here r is related to t by e H(,t = —Hot. Then t — > —00 corresponds to r 
As an example of f(r) in (|73p . we may consider 



f(r) 



HIt 2 (1 + H]t 2 ) ' 



(76) 



(77) 



-00 and t — > +00 corresponds to r — > 0. 



(78) 
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where Hl and Hi are constants. At early times in the history of the universe r — > oo (corresponding to t — > — oo), 
/(t) behaves as 

/(T) lp • (79) 

Then the Hubble rate is given by a constant Hi, and therefore the universe is asymptotically de Sitter space, corre- 
sponding to inflation. On the other hand, at late times r — > (corresponding to t — > +00), /(r) behaves as 

/w - #p ■ ( 8 °) 

Then the Hubble rate is again a constant Hl, which may correspond to the late time acceleration of the universe. 
This proves that our reconstruction program can be applied directly to the non-minimally coupled scalar theory. 



V. DE SITTER SPACE IN SCALAR-TENSOR THEORY 



When studying scalar field cosmology in a spatially flat FRW universe from the dynamical systems point of view, 
it is often convenient to redefine the scalar field <p used in the previous sections in such a way that its kinetic energy 
density has a canonical form (apart from the sign). Instead of the field <fr appearing in Eq. (|47j) . one can use 



(81) 



in terms of which the action (l47l) becomes 



S = 



d 4 x. 



-9 



(82) 



where e = sign(JT) and U(a) = U [</>(<r)] (this redefinition, however, can not change the sign of the kinetic energy 
of the phantom field to make it positive - for this purpose one would have to make the scalar field purely imaginary 
through a sort of Wick rotation). 

In discussions of the phase space of spatially flat scalar field cosmology, the Hubble parameter H and the scalar 
field a (or (f>) constitute a natural choice of dynamical variables. The phase space is a two-dimensional curved surface 
embedded in a three-dimensional space and this suffices to guarantee the absence of chaos in the dynamics [l3|. 
Moreover, the only fixed points of the dynamical system are de Sitter spaces with constant scalar field (Hq, <jq). For 

the solution described by Eqs. (|54|) and ([55]) . the scalar field redefinition (|8 1 [) yields a =const.= <jq and (^Ho,a<jj 

is a fixed point of the system. The fixed point nature of a particular solution does depend on the specific choice of 
dynamical variables: for example, the solution (T54")l and {55]) is a fixed point with the choice (H, a) but not with the 
choice (H,4>). While the solution is the same, it is convenient to study the dynamics using a instead of cf>. This is 
particularly important for detailed calculations of the stability of de Sitter space using gauge-invariant variables (see 
below), as these calculations greatly simplify in a de Sitter background. 

One may wonder whether the fixed point (or even the attractor) nature of de Sitter spaces is lost when performing 
the conformal transformation to the Einstein frame used in the previous sections in order to find exact solutions. This 

is not the case, but a little care is needed because, in general, acceleration of the universe a > in the Jordan frame 

.2 ~ 

does not imply cosmic acceleration > in the Einstein frame, and vice- versa since 



(83) 



(here an overdot denotes differentiation with respect to the Jordan frame comoving time t). However, a de Sitter space 
in the Jordan frame is mapped into a de Sitter space in the Einstein frame (and vice- versa). A general conformal 
transformation of the metric g^ LV — » g^ v — £! 2 g M „ , where the conformal factor fl depends on the scalar field present 
in the theory (for example, as in Eq. (|45| ). yields the rescaled FRW line element 



ds 2 = n 2 ds 2 = -di 2 + a 2 (t) (dx 2 + d y 2 + dz 2 ) 



(84) 
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with dt = Sldt and 5 = £1 a. Therefore, the relation between the Hubble parameters of the two conformal frames is 























a dt 




" n J 







(85) 



Since f2 = £l((r), it is clear that a Jordan frame fixed point yH, a J — (0,0) is mapped into an Einstein frame fixed 

point ^iJ, cr^ = (0,0). Moreover, a small perturbation of this fixed point in the Jordan frame corresponds to a 

small perturbation in the Einstein frame, and stability of the Jordan frame fixed point corresponds to stability of the 
corresponding Einstein frame stationary point. In fact, assume that 

H(t) = Ho + SH(t) , (86) 



a(t) = <jq + Sa(t) , 



(87) 



in the Jordan frame, where \5H(t)/Ho\ and \8a(t)/ao\ are small perturbations which, for simplicity, are taken here to 
be homogeneous. Then, in the Einstein frame, H = Hq + SH with 



SH 



1 



SH +7T ( S& - H oH 
"0 



(where a prime denotes differentiation with respect to a and the zero subscript denotes quantities evaluated in the 
background de Sitter space). The Einstein frame perturbation stays small if the Jordan frame perturbations are 
small, hence an attractor in the Jordan frame corresponds to an Einstein frame attractor (the converse is not always 
true, see [l4| for a counterexample). This property is crucial when using conformal transformations to study slow-roll 
inflation which describes dynamics around a de Sitter attractor. It is the presence of a de Sitter attractor that justifies 
the use of the slow-roll approximation H(t) — Hq + 8H(t) in inflation, and the fact that the presence and attractor 
nature of de Sitter space (subject to certain conditions) are guaranteed also in the Einstein frame ultimately justifies 
the use of conformal techniques in the study of slow-roll inflation. To summarize, in general, cosmic acceleration in 
one conformal frame does not correspond to acceleration in the conformally related frame(see the related discussion 
in [HI]), however slow-roll inflation in the Jordan frame corresponds to slow-roll inflation in the Einstein frame. This 
is relevant also for late time de Sitter-like expansion in a universe dominated by quintessence. 

The previous discussion is restricted to homogeneous (space-independent) perturbations of de Sitter space, but it 
can be extended to more general (space-dependent) inhomogeneous perturbations. The latter are more problematic 
because of the notorious gauge-dependence problems associated with them, and they are best described using gauge- 
independent methods. A linear stability condition of de Sitter space against inhomogeneous perturbations was derived 
in [161 ] . This is a necessary condition for de Sitter space to be an attractor in the (H, a) phase space, however it is 
not sufficient because it only ensures stability to first order in the gauge-invariant variables (stability to higher, or to 
all orders, can usually be established only numerically and with the restriction to homogeneous perturbations, see, 
e.g., O). 

For a theory described by the action 



S : 



the gauge-invariant linear stability condition is [16 



JE I 
d> 2 lo 



d-*V 



b 2 io 



+ 



Fo 



u H3 



loo Fo 



< 



(89) 



(90) 



where F = ^ and (p,pji= -§^§^- For the action (|44|) . this condition becomes 



/o'(i + /o) + (/oy 



2 K ^o (i + fo) + 3 (/ r 



< 



(91) 



where f = f{(fio)>fb = f'{4>o)> e ^ c - an d assuming that 1 + f(c/>o) > in order to guarantee a positive effective 
gravitational coupling. For an "ordinary" phantom field in general relativity it is <p(4>,R) = R, uj — —1, and the 
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condition (|9T|) reduces to Vq < 0, i.e., the de Sitter space is stable if the potential has a maximum to which the 
phantom field can climb and settle in during the dynamical evolution. This behaviour, which is opposite to that of 
an ordinary scalar field, is due to the negative sign of the kinetic energy of the phantom 0]. 
The FRW equations give two conditions for the existence of de Sitter fixed points: 



Ro || | = 2^o - V , (92) 



= Vo • (93) 



For the action (|44|) these become 



VS = \±fv , (95) 

J- + JO 

where Rq = 127Jq for de Sitter space. Upon use of Eq. (|M|) , the stability condition (f9"Tj) in this theory becomes 

(U / '-2y ")(l + / ) + U (^) 2 
2^u, (l + / ) + 3(/ ) 2 " 

These stability conditions are important in the study of the future evolution of the late-time accelerated era. 



VI. LATE TIME ACCELERATION AND INFLATION WITH SEVERAL SCALAR FIELDS 



In this section we begin by considering a model with two scalar fields minimally coupled to gravity (see pi 111 ) . 
Such models are used, for example, in reheating scenarios after inflation. 

An additional degree of freedom appears in this case, so that for a given solution we may choose different conditions 
on the scalar fields, as shown below. It is possible to restrict these conditions by studying the perturbative regime for 
each solution. The action we consider is 



S = 



A _ oj{4>)d^m - a( x )d^ l x - V(4>, X) 



(97) 



where u>(4>) and <j(x) are the kinetic terms, which depend on the fields <j> and x> respectively. We again assume a flat 
FRW metric. The Friedmann equations are written as 



H 2 



-2 r 



1 



2 + \<r(x)x 2 + V(<f>,x) 



H = 



K 

Y 



uj{(t>)4> 2 + <j(x)x 2 



By means of a convenient transformation, we can always redefine the scalar fields so that we can write <j) 
The scalar field equations are given by 

" + + 3Hu(<f>)j> + = , a{ X )x + \°'{X)X 2 + ZHa{ X )x + = 



(98) 

x = t. 

(99) 



Then, for a given solution H{t) = f(t), and combining the first Friedmann equation with each scalar field equation, 
respectively, we find 



,q = 2 df{4>, x) , )= 2 0/(0, x ) 

k 2 dtp k 2 dx 

where the function /(</>, x) carries down to f(t,t) = f(t), and is defined as 



KM 



K 

Y 



w(^)d0+ / <r(x)dx 



(100) 



(101) 
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The scalar potential can be expressed as 



v{4>, x ) 



and the second Friedmann equation reads 



3/(0, x) 



2 , W, X ) , 5/(0, X ) 



(102) 



,f'{t) = oj{t)+a{t). 



(103) 



Then, the kinetic functions may be chosen to be 

"(0) = -4 TO) + am , <K0 = > 

where 5 is an arbitrary function. Hence, the scalar field potential is finally obtained as 

1 



v(4>, x ) = - [3/(0, xf + /'(0) + <?(0) - g(x)j 



(104) 



(105) 



A. Example 1 



We can consider again the solution (fT7|) 



/(*) 



Hq Hi 



t.-t t 2 



(106) 



This solution, as already seen in Sec. II, reproduces unified inflation and late time acceleration in a scalar field model 
with matter given by action JT]). We may now understand this solution as derived from the two-scalar field model 
(|97| . where a degree of freedom is added so that we can choose various types of scalar kinetic and potential terms, as 
shown below for the solution p06p . Then, from Eqs. Q104p and (|106[) . the kinetic terms follow: 



Ho 



(ts - 0) 2 



2 —+a{4>) 



> U (4>) = -ndiX) 



The function f(4>,x) i s 



M,x) 



"" +§ + / #5(0)- / d X g(x), 



ts - 



(107) 



(108) 



while the scalar potential is 



3/(0, xY 



Hn 



2H 



(ts - 0) 2 ; 



1 + 9(4) - 9<X) 



(109) 



This scalar potential leads to the cosmological solution (fT7]) . 

It is possible to further restrict g(t) by studying the stability of the system considered. To this end, we define the 
functions 



X x = x, Y 



/(0,X) 
H 



Then, the Friedmann and scalar field equations can be written as 

dX^ _ 1 c/(0) lv2 



dN 
dX a 

dN 2H <j( X ) 
dY k 2 



2H^V^-l)-3(^-y). 
1 a ' ix) (Xl-1)-3(X X -Y) , 



(110) 



(111) 



dN 2H 2 



■ufflX+prXi, - 1) + a( X )X x (YX x - 1)] , 



(112) 
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where 



377 = Tilt- At = X x = Y = 1, we consider the perturbations 



then 




Xa = 1 + <DQ 



a x = i + £a x , y = i + sy, 



= M I JA X 

jy 



M = 









*'(*) _ 

Ha( X ) s 

,..2 q-(x) 

K 2ff 2 



(113) 



(114) 



,2 "(0)+o-(x) 
2ff 2 



The eigenvalue equation is given by 



3 + A 



3« 2 w(0) / <7'( X ) 



+3 + A^) f-^( w (0) + (T (0))-A 



+ 3 + A + 



V2# 2 



3 K 2 a( X ) / 



2H 2 \Ha(x) J 2H 2 \Hu(4>) 

To avoid divergences in the eigenvalues, we choose the kinetic functions to satisfy 

w(4>) + , a( X ) ± , 



3 + A = 



(115) 



(116) 



hence, the eigenvalues in Eq. pi5[) are finite. Summing up, under these conditions, the solution (|106[) has no infinite 
instability when the transition from the non-phantom to the phantom phase occurs. 

As an example, we may choose g(t) — a/t 3 , where a is a constant that satisfies a > 2H\. Then, the /(</>, x) function 
(|108|) is given by 



f(<i>,x) 



Ha 



t s - 4> 



{a-2H 1 ) 
24 2 



As a result, the kinetic terms fllOTf) are expressed as 



and the potential reads 



v(<p,x) = i 



Ho 



2H 1 a 



% ~ 0) 2 3 



3/(</>,x) 5 



a 
2^ 



2 a 

cr{X) = — — 
K z X 



Hn 



a-2H x 



a 



(117) 



(118) 



(119) 



This potential reproduces the solution above, which unifies inflation and late time acceleration in the context of 
scalar-tensor theories, involving two scalar fields. Notice that the extra degree of freedom gives the possibility to 
select a different kinetic and scalar potential in such a manner that we get the same solution. 

In the case in which the condition (|116p is not imposed, the kinetic terms (| 104[) may have zeros for < t < t s , so 
that the perturbation analysis performed above ceases to be valid because some of the eigenvalues could diverge. 



B. General case: n scalar fields 



As a generalization of the action (|97p . we now consider the corresponding one for n scalar fields, 



S = J d A x^—g 
The associated Fricdmann equations are 



f r i 1 - 

J d 4 Xy^ ^-j^WW^-^lfe'" An) 



H z = — 



n 1 



UJi{<t>i)(t)i + V((f>i, ■■■ ,4>n) 



2 



(120) 



(121) 
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We can proceed analogously to the case of two scalar fields so that the kinetic terms are written as 

n „ 

5>i(t)=--5/'(t), (122) 



K 
i=l 



hence, 



Ui{(Ri)=--2 7J i y{.9u--'i9n) = -2 



3/(0i, ■■■ ,<j>n) +}^ TT 



i=l 



(123) 



where /(t, i, • • • , t) = f(t). Then the following solution is found 

9i = t , H(t) = f(t) . (124) 
From (|122p we can choose, as done above, the kinetic terms to be 

2 2 2 

Wl(</>l) = o I/'(<M + 52(0l) H 1" ffn(0l)] , ^2(^2) = ^32(02) , ' ' ' , W„(0„) = -zgnipn) ■ (125) 

ft ft ft Z 

Then, there are n — 1 arbitrary functions that reproduce the solution (|124|) so reconstruction may be successfully 
done. They could be chosen so that dark matter is also represented by some of the scalar fields appearing in the 
action ([120]) . 

VII. RECONSTRUCTION OF INFLATION AND COSMIC ACCELERATION FROM TWO-SCALAR 

THEORY 

In the present section, inflation and cosmic acceleration are reconstructed separately, by means of a two scalar field 
model that reproduces some of the cosmological constraints at each epoch. We explore an inflationary model in which 
the scalar potential, given for a pair of scalar fields, exhibits an extra degree of freedom and can be chosen in such 
way that slow-roll conditions are satisfied. Also, the cosmic acceleration is reproduced with a pair of scalar fields plus 
an ordinary matter term, in which the values of the observed cosmological density parameter (Ode — 0.7) and of the 
EoS parameter (wde — — 1) are reproduced in a quite natural way. For the case of a single scalar, the reconstruction 
for similarly distant epochs was given in Ref. 19] (for other reconstruction versions, see also [20I]). 



A. Inflation 



In the previous sections, models describing inflation and late-time accelerated expansion have been constructed by 
using certain convenient scalar-tensor theories. In this section, we present an inflationary model with two scalar fields, 
which can be constructed in such a way that the inflationary conditions are carefully accounted for. For this purpose, 
we use some of the techniques given in the previous section. The action during the inflationary epoch is written as 



S 



-fJ 



R 
2^2 



1 



1 



(126) 



We will show that a general solution can be constructed, where the scalar field potential is not completely specified 
because of the extra degree of freedom represented by the second scalar field added, in a way similar to the situation 
occurring in the previous section. Considering a spatially flat FRW metric, the Fricdmann and scalar field equations 
are obtained by using the Einstein equations and varying the action with respect to both scalar fields: 



H 2 = — 
3 



2^ + 2 X 



• + 3H(j> + 



, X + ZH X - 



dx 



. 



(127) 



We assume that the slow-roll conditions are satisfied, that is, <f> <C 3H<j> (jl "C 3Hx), and <j> 2 -C V(cf), x) (x 2 ^ ^('Aj X)) > 
in order for inflation to occur. Then, Eqs. (|127[) take the form 



H- 



-V(4>, X ), 3iJ0 + 



dV(4>, X ) 



0, 3H X 



dV(y,x) 
dx 







(128) 
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and the slow-roll conditions read 



1 ViV> 1 1 JViiVv , s 

• <1, 77- — 77 <1. (129) 



3k 2 V 2 ' 3k 2 V 

Here, Vj denotes the partial derivative of V with respect to one of the scalar fields (i = <fi,x)- As done previously, a 
scalar potential V((f>,x) can be constructed, although in this case the conditions for inflation need to be taken into 
account. From (|128|) , the potential is given by 

V{<t>, X ) = \h 2 {^ X ) ■ (130) 

We can choose this potential such that 

V(^x) = ^[f 2 (<P,x)+9i(^-92(x)] ■ (131) 

The three components /, g%, and gi are arbitrary functions, and g(N) = gi(N) = g%{N) where, for convenience, we 
use the number of e-folds N = In instead of the cosmic time, and a* denotes the initial value of the scale factor 
before inflation. Then, the following solution is found: 

H(N) = f(N) . (132) 

Hence, Eqs. (|128|) may be expressed as a set of differential equations with the number of e-folds as independent 
variable, 

To illustrate this construction, let us use a simple example. The following scalar potential, as a function of the 
number of e-folds N, is considered: 

V{^ X ) = \[HlN 2a ] , (134) 

where a and Hq are free parameters. By specifying the arbitrary function g(N), one can find a solution for the scalar 
fields. Let us choose, for the sake of simplicity, g(N) = goN 2a , where go is a constant, and f((f>,x) = /(</>) (i-e., as a 
function of the scalar field only). Then, using Eqs. (|133p . the solutions for the scalar fields are found to be 

0(iV) = <f> - -4rJMH§ +90) hxN , x(N)= X o-^rV2g^N , (135) 
kHo v kHo 

and the scalar potential can be written as 



,„ , v ( V^Hq \ fn 2 H 2 ( X0 -xf 

(H + g Q ) exp k (<pa - <p) - go' 



We are now able to impose the slow-roll conditions by evaluating the slow-roll parameters 

1 VjV* _ 2a ({Ho+go? , 4 ff0 , 



(136) 



3k 2 V 2 3H 2 \ H 2 +g N 



1 V^W 71 2 2 (Hp+go) 2 16,g (4 Q -l) 2 1 

~ > a \ 7H2 , „ N 2 + — 77273 — nn^ 1 - ( 137 ) 



3k 2 y 3 V (H 2 + go) 2 4g 2 a 2 N 2 

Hence, we may choose conveniently the free parameters so that the slow-roll conditions are satisfied, and therefore, 
inflation takes place. From these expressions we see that the desired conditions will be obtained, in particular, when 
a is sufficiently small and/or Ho and N are large enough. All these regimes help to fulfill the slow-roll conditions, in 
a quite natural way. 
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B. Cosmic acceleration with a pair of scalar fields 



It is quite reasonable, and rather aesthetic, to think that the cosmic acceleration could be driven by the same 
mechanism as inflation. To this purpose, we apply the same model with two scalar fields, with the aim of reproducing 
late-time acceleration in a universe filled with a fluid with EoS p rn — w rn p m . The free parameters given by the model 
could be adjusted to fit the observational data, as shown below. We begin with the action representing this model, 



-9 



R 
2^2 



By assuming a spatially flat FRW metric, one obtains the Friedmann equations 



H 2 = — 



\(4>) 2 + \(x) 2 + v(^x)+p« 



H 



\<t> 2 + \x 2 + v(4>, x ) + Pn 



(138) 



(139) 



Variation of the action (|138p yields the scalar field equations 

4> + 3H4> + = , X + "SHx + V x = . (140) 
This set of independent equations may be supplemented by a fifth one, the conservation of matter-energy density, p m , 



3Hp m (l + w m ) = 



As done in [19j , we perform the substitutions 



3H 2 /k 2 



sf 



H<j>) 2 + h(x) 2 + v(^x) = H_ 

3H 2 /k 2 ' 6 H 2 



IJ sf 



p sf h^) 2 +m 2 ~y^x) 

Psf \$f + \{x)* + V{hx) 



(141) 



(142) 



For convenience, we consider both scalar fields together under the subscript sf, so that the corresponding density 
parameter is ft s f = il c f, + il x . Then, using the transformations Q142|) , the Friedmann equations and the scalar equation 
read 



n s/ = i, 2e + 3(i + w sf )n sf + 3{i + w m )n m = , n' sf + 2en sf + sn sf (i 



= 



(143) 



Here, the prime denotes differentiation with respect to the number of e-folds N = In 



We can now combine 



Eqs. p42p to write the EoS parameter for the scalar fields as a function of tt s f and of the time derivative of the scalar 
fields, namely 



K 2 W>' 2 + X ,2 )-3fi s/ 



sf 



(144) 



Hence, it is possible to find an analytic solution for the equations (|143p . for a given evolution of the scalar fields, as 
will be seen below. Before doing that, it is useful to write the effective EoS parameter, which is given by 



i 2e 

l"cff = -1 - y , Pcff = Pm + Psf , PcS = Pm + Psf , 



(145) 



and the deceleration parameter 



q = 



iH 2 



= -1 - e 



(146) 



As usual, for q > the universe is in a decelerated phase, while q < denotes an accelerated epoch, such that 
for w c h < —1/3 the expansion is accelerated. To solve the equations, we consider a universe which, at present, is 
filled with a pressureless component (w m = 0) representing ordinary matter, and two scalar fields which represent a 
dynamical dark energy and a dark matter species. To show this, we make the following assumption on the evolution 
of the scalar fields, which are given as functions of N: 



(147) 
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Then, Eqs. (] 143[) can be solved, and the scalar field density parameter takes the form 

Vs f = ^ + n x = i- 1 ^— , (W8) 

where k is an integration constant and A = 3 — (a 2 + /3 2 ). It is possible to introduce an arbitrary function g(N), to 
express the energy density parameter for each scalar field in the following way: 



il = ^- ke XN , V x =g(N). (140, 



The function g(N) may be chosen in such a way that the scalar field \ represents a cold dark matter contribution at 
present (w x ~ 0), and the scalar field ef> represents the dark energy responsible for the accelerated expansion of our 
universe. On the other hand, using Eqs. (|143|l . e = Mr is obtained as 



kX (ke XN + a 2 /3 2 



(150) 



"kg = - 1 - o £ = - r/-2 , 2N — AJV Z Ul (h^N I o \ ' ( 151 ) 



[(a 2 + /3 2 )e- AAr + 3k] (ke XN + 3) J 
Then, it is possible to calculate the effective parameter of EoS given by Eq. (| 1 33[) 

2 fcA (/fce AAr + a 2 (3 

3 £ ~ ~ [(a 2 + p 2 )e~ XN + 3k] (ke XN + 3) 



We have four free parameters (N, k, a, (3) that may be adjusted to fit the constraints derived from observations. 
With this purpose, we use the observational input fi TO ~ 0.03, referred to baryonic matter, and normalize the number 
of e- folds N, taking N = at present, then the integration constant k may be written as a function of a and f3 as 

n m (N - 0) - 0.03 -> k = 2 " 01 ~ ( o3 + /?2) ■ (152) 

The free parameter f3 may be fixed in such a way that the scalar field \ represents cold dark matter at present, i.e., 
w x ~ and f2 x ~ 0.27, and its EoS parameter is written as 

_ n 2 X ' 2 - 3il x _ P 2 - 3g{N) 
WX ~ 30 x " 3g(N) ■ (153) 

For convenience, we choose g(N) = ^-e~ N , then the energy density and EoS parameter are given by 

92 



1 — P~ iv If 

™x = —^: «x = ye- w . (l.vl :) 



Hence, at present (N — 0), the expressions (|154[) can be compared with the observational values and the (3 parameter 
is given by 

w x (N = 0) = 0, n x {N = 0) ~ 0.27 -► 1 = 0.81 . (155) 
Finally, the energy density of (j> expressed by Eq. (|149[) takes the form 



2 



The value for a could be taken so that ~ 0.7 and w<p ~ — 1 at present. Hence, it has been shown that cosmic 
acceleration can be reproduced with a pair of scalar fields, where due to the presence of the extra scalar that can be 
identified with the dark matter component. 

It is interesting to point out that one can unify these realistic descriptions of the inflationary and late-time accelera- 
tion eras within a single theory. However, the corresponding potential looks quite complicated. The easiest way would 
be to use step (^-function) potentials in order to unify the whole description in the easiest way (as was pioneered in 

We may also construct a model unifying early universe inflation and the present accelerated expansion era. To 
this end we can choose f{<f>) in (f30|) . which gives the Hubble parameter (|3T|) . Then, using (|104[) . one can define oj(<j>) 
and cr(x) with the help of an arbitrary function g. After defining then f(4>,x) with Eq. (jlOip . we can construct the 
potential V(<p, x) using Eq. (|102p . Finally, we obtain the two scalar- tensor theory (f9"T)l reproducing the Hubble rate 
(|3"Tj) . which describes both inflation and the accelerated expansion. 
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VIII. CONCLUSIONS 



Modeling both the early inflation and late-time acceleration epochs within the context of a single field theory has, 
undoubtedly, much aesthetic appeal and seems a worthy goal, which we have attempted here. To summarize, we have 
developed, step by step, the reconstruction program for the expansion history of the universe, by using a single or 
multiple (canonical and/or phantom) scalar fields. Already in the case of a single scalar, we have presented a number 
of examples which prove that it is actually possible to unify early-time inflation (at very high redshift) with late-time 
acceleration (at low redshift). The reconstruction technique has then been generalized to the case of a scalar field 
non-minimally coupled to the Ricci curvature, and to non-minimal (Brans-Dicke-type) scalars. Again, various explicit 
examples of unification of early-time inflation and late-time acceleration have been presented in those formulations. 
Due to the special role of de Sitter space, which often appears as an attractor in the inflationary epoch, as well as in 
the present cosmic acceleration era, special attention was paid to this specific space-time. Conformal transformations 
to the Einstein frame and stability conditions for the de Sitter space were discussed. 

Moreover, the case of several minimally coupled scalar fields has been considered for the description of the realistic 
evolution of the Hubble parameter, and we have shown that it is qualitatively easier to achieve a realistic unification 
of late and early epochs in a model of this kind, in such a way as to satisfy the cosmological bounds coming from the 
observational data. This is due to the arbitrariness in the choice of the scalar potential and the scalar kinetic factor 
in the description of a universe with a given scale factor a(t). Using the freedom of choosing these scalar functions, 
one can constrain the theory in an observationally acceptable manner. Specifically, slow-roll conditions and stability 
conditions may be satisfied in different ways for different scalar functions, while the scale factor remains the same. 
This can be used also to obtain the correct structure of perturbations, etc. 

As a mater of fact, many questions remain to be discussed in greater detail, a more realistic matter content should 
be taken into account, and the universe expansion history be described in a more precise and detailed manner. After 
all, we live in an era of increasingly more precise cosmological tests. Anyhow, the unified effective description of 
the cosmic expansion history presented here seems quite promising. Using it in more realistic contexts — in which, of 
course, technical details become necessarily more complicated — appears to be quite possible. 
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